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ABSTRACT 

/  ^ 

r  Heat  release  analysis  can  be  computed  from  the  measured  cylinder  pressure  at 
known  cylinder  volumes  through  the  use  of  the  First  Law  of  Thermodynamics  and 
the  ideal  gas  law.  After  correcting  for  crevice  effects  and  heat  transfer  to 
the  combustion  chamber  walls,  the  combustion  heat  release  profile  can  be 
determined.  From  the  combustion  heat  release  profile,  the  burn  rate  and  the 
mass  fraction  burned  can  be  calculated  for  studies  of  knock  and  turbulence. 

The  Ricardo  Hydra  MK  III  engine  is  instrumented,  the  operating  properties  are 
checked  and  the  cylinder  pressure  is  measured  for  fired  and  motored  operating 
conditions  at  six  different  load  and  speed  setpoints.  The  generated  pressure 
data  is  then  used  in  the  Sloan  Automotive  Lab  heat  release  model  to  determine 
the  validity  of  the  model  for  real  engine  data.  The  model  is  checked  for 
scaling  by  load  and  speed.  The  heat  release  was  successfully  modeled  using 
one  set  of  model  parameters  for  this  engine  operating  at  the  different  load3 
and  speeds. 
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CHAPTER  ONE 


INTRODUCTION 


Many  studies  have  been  made  using  internal  combustion  engine  pressure  data 
to  determine  engine  preformance.  The  thermal  energy  of  the  working  fluids 
inside  the  engine  can  be  accounted  for  through  use  of  the  cylinder  pressure 
data.  Heat  release  analysis  can  help  determine  the  energy  release  during 
combustion  and  the  energy  loss  from  heat  or  mass  transfer  during  an  engine 
cycle.  Few  studies  have  been  made  to  determine  how  well  the  heat  release 
model  will  scale  with  engine  load  and  engine  speed.  A  study  of  this  nature 
would  assist  in  determining  the  validity  of  a  heat  release  model. 

Heat  release  analysis  is  computed  from  measured  cylinder  pressure  at  known 
cylinder  volumes  and  through  use  of  the  First  Law  of  Thermodynamics  and  the 
ideal  gas  law.  This  allows  for  determination  of  the  apparent  heat  release 
Inside  the  cylinder.  By  adding  the  effects  of  the  heat  transfer  to  the  walls 
and  the  affects  of  the  cylinder  crevices,  the  heat  release  due  to  combustion 
can  be  determined.  The  value  of  the  total  heat  released  during  combustion  can 
be  compared  to  the  chemical  energy  introduced  for  that  cycle  to  determine  the 
validity  of  the  model. 

This  work  deals  in  three  main  areas.  First,  the  engine  setup, 
instrumentation  and  verification  checks  for  the  Ricardo  HYDRA  MK  III  engine 
from  which  the  pressure  data  was  taken  and  reviewed.  Second,  the  data 
aquisition  procedure  is  described.  Third,  the  heat  release  analysis  and 


results  are  discussed 


CHAPTER  TOO 


THEORETICAL  BACKROUND 


2.1  FIRST  LAW  ANALYSIS 

Analysis  of  the  cylinder  contents  of  an  internal  combustion  engine  usually 
starts  with  the  first  law  for  an  open  system. 

«U  -  6Q  -  6Vf  ♦  Jln1h1 

where  U  is  the  internal  energy  of  the  gases  Inside  the  cylinder,  Q  is  the  heat 
transfer  across  the  boundaries,  W  is  the  work  done  by  the  system  on  the 
control  volume  and  is  the  energy  flow  due  to  mass  flow  across  the  system 
boundary. 

For  the  internal  combustion  engine  system,  U  represents  the  sensible 
internal  energy  which  is  a  function  of  temperature  and  can  be  written  as: 

U  -  mcyT 

where  m  is  the  mass  in  the  cylinder,  cy  is  the  constant  volume  specific  heat 
and  T  is  the  average  cylinder  temperature.  Using  the  ideal  gas  law  pV  -  mRT, 

U  can  be  rewritten  as: 


The  piston  work,6W,  can  be  written  as  pdV. 

The  heat  transfer  Q  can  be  divided  into  two  parts: 

Q  ,  the  heat  addition  due  to  combustion  and 
c 

Qy,  the  heat  transfer  to  the  combustion  chamber  walls. 

Using  the  additional  relations:  c-c»R,Y-c/c  ,  h  -  c  T  and 

p  v  p  v  p 

rearranging,  the  first  law  equation  becomes: 
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SQc "  5Qw +  y^t  Vdp  +  Vdp  +  cpTdm 

Considering  a  closed  system  (figure  2-1 )  from  intake  valve  closing  to 
exhaust  valve  opening,  eliminates  the  changing  mass  term.  The  first  law 
equation  becomes: 

«QC  "  5Qw  +  TT  VdP  +  y=T  VdP 

The  net  heat  release,  equal  to  the  heat  added  by  combustion  minus  the  heat 
transferred  to  the  walls,  can  be  calculated  from  the  measured  pressures,  the 
known  cylinder  volune,  and  the  known  rate  of  change  of  the  cylinder  volume. 


2.2  CREVICE  'AFFECT 

For  the  first  law  derivation  of  the  heat  release  equation,  the  system  was 
considered  to  be  a  closed  system  and  the  changing  mass  term  was  dropped. 
However,  there  are  spaces  with  low  volume  to  surface  area  ratios  where  the 
wall  temperature  has  a  greater  efect.  Usually  called  crevices,  these  spaces 
include  gaps  for  the  head  gasket,  for  the  spark  plug  threads,  and  the  volume 
between  the  piston  and  the  cylinder  wall  down  to  the  piston  rings  (which  is 
typically  the  largest  crevice  volume).  These  regions  can  be  grouped  together 
and  modeled  as  a  single  crevice  control  volume  (figure  2-1).  The  crevice 
volume  is  typically  1  to  2%  of  the  clearance  volume,  but  due  to  the  cooling 
effect  of  the  large  surface  area,  the  crevice  may  contain  6-10%  of  the 
cylinder  mass. 

The  temperature  of  the  mixture  in  the  crevice  can  be  modeled  as  equal  to 
the  temperature  of  the  cylinder  gases  for  mass  flow  into  the  crevice  but  equal 
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Figure  2-1.  Cylinder  and  Crevice  Control  Volumes 


to  the  wall  temperature  for  mass  flow  out  of  the  crevice.  This  difference  in 
temperature  of  the  mass  flew  represents  a  loss  of  energy  from  the  cylinder 


control  volume.  To  account  for  this  energy  change,  the  mass  effect  term  for 

the  crevice  must  be  included  in  the  heat  release  equation.  This  term: 

c  Tdm 
P 

can  also  be  changed  by  using  the  relation  pV  -  mRT  as  follows: 

Vdp  * 

cpT(dB)*CpT[lf]-^|  V* 

cr  cr 

* 

where  T  -  crevice  wall  temperature  for  d  >  0  and 

m  — 

* 

T  ■  cylinder  gas  temperature  for  d  <  0 

m 

cr  »  indicates  crevice 

With  this  term,  the  heat  release  equation  becomes: 


j  V<p 

cr 


2.3  HEAT  TRANSFER 

An  estimate  for  the  heat  transfer  to  the  walls  must  be  made  if  the 
combustion  heat  release  is  to  be  separated  frcm  the  net  heat  release.  In 
general,  heat  transfer  from  a  hot  gas  to  the  cool  combustion  chamber  surface 
is  written  as: 

Q  -  hAAT 

where  h  is  the  heat  transfer  coefficient,  A  is  the  combustion  chamber  surface 
area  and  AT  is  the  temperature  difference  between  cylinder  gas  and  the 
cylinder  wall.  Woschnl  (1)  developed  an  experimental  correlation  for  heat 
transfer  based  on  the  Nusselt-Reynolds  relationship: 


The  characteristic  length  used  is  the  bore  while  the  characteristic  velocity- 
in  the  Reynolds  Number,  Re,  is  proportional  to  the  mean  piston  speed  with  a 
factor  added  for  the  combustion  expansion  velocity.  This  factor  is  used 
during  the  combustion  and  expansion  phases  of  the  engine  cycle. 

Woschni's  correlation  for  the  heat  transfer  coefficient  is: 

h(w/m2k)  -  ^BGnf0,2  p(atm)0-8  TOC)-0,53  W(m/s)0,8 

For  this  equation,  c1  is  a  constant,  B  is  the  cylinder  bore,  p  and  T  are  the 
average  cylinder  pressure  and  temperature  and  W  is  the  characteristic 
velocity.  The  value  for  the  characteristic  velocity  is  calculated  from: 


where  Sp  is  the  mean  piston  speed,  c2  is  a  constant,  Vp  is  the  cylinder 
volume,  Pf  and  Pffl  are  the  fired  and  the  motored  pressures  for  that  position, 
and  V^.P^.and  T^  are  the  cylinder  conditions  at  the  intake  valve  closing. 
Woschni  empirically  found  c1  to  equal  110  and  c2  to  equal  3.24  X  10  3.  To 
achieve  consistancy  with  S.I.  units,  the  value  for  c1  was  changed  to  3.26  in 
the  computer  program. 

For  the  engine  anaylsis,  the  computer  program  used  to  calculate  the  heat 
release  was  written  to  multiply  Woschni's  c1  (-3.26)  and  c2  (-3.24  X  10~3)  by 
an  input  value  of  C,  and  0,.  The  input  value  of  C,  and  C2  are  the  variables 
that  scale  the  heat  release  model  to  the  Ricardo  engine  data. 


Woschni ' s  heat  transfer  correlation  can  also  be  modified  by  changing  the 
exponents  which  scales  the  heat  transfer  coefficient  with  engine  speed. 
Woschni' s  formula  can  be  written  as: 

h-o,  pnT(0-re-,-62n,»n  (1) 

where  n  is  the  exponent.  With  n  -  0.8,  the  original  equation  is  recovered. 

2.4  DETERMINATION  OF  GAMMA  (Y) 

Y  (-Cp/cv)  is  an  important  parameter  in  the  heat  release  model.  It  is  a 
parameter  which  varies  as  a  function  of  both  pressure  and  temperature,  and 
with  the  proportions  of  unbumed  and  burned  gases.  To  use  y  in  the  heat 
release  model,  it  is  convenient  to  regard  Y  as  a  single  function  of 
temperature,  in  this  case,  a  linear  function  of  temperature. * 

To  obtain  Y  as  a  function  of  temperature,  cylinder  pressure  and 
temperature  were  calculated  from  the  experimental  data  at  various  points  and 
plotted  on  the  chart  of  Y  versus  temperature  as  given  in  Mansouri  and 
Heywood(2).  A  calculation  of  Y  for  the  unbumed  mixture  was  made  for  the 
temperature  from  500°K  to  1500°K.  These  points  were  also  plotted  and  a  least 
squares  linear  curve  fit  was  drawn  through  both  sets  of  points  (figure  2-2). 
The  equation  for  this  line  was  incorporated  into  the  heat  release  model. 

The  mass  in  the  cylinder  was  required  in  order  to  calculate  the 
temperature  in  the  cylinder  using  the  ideal  gas  law.  The  mass  in  the 
cylinder: 


as 


Figure  2-2.  y  vs.  T  for  Indolene  at  $  -  1.0 
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mass  -  m  ♦  m,  + 
a  r 


can  be  written  as: 


m 


,4  r  *3  [f l  120  .  «r  tf )  12°  1 

1  n  RPM  RPM 


where  Is  the  residual  fraction,  is  the  mass  of  air,  and  mf  is  the  mass 
of  fuel.  Using  the  measured  pressure  and  the  known  volume: 

t  .  p  *  TOL 
m  x  R 

The  plotted  values  for  Y  unburned  were  calculated  using  a  mass-fraction 
weighted  sun  of  Y  air  and  Y  fuel  at  $  •  1 .0.  Y  fuel  was  determined  by 
interpolating  Indolene's  average  molecular  weight  between  the  value  of 
and  CgH^.  The  mass  weighted  value  of  Y  for  the  fuel  air  mixture  was 
subtracted  for  Y  air  and  the  difference  was  plotted  on  the  Y  versus 
temperature  chart.  The  result  of  the  least  squares  curve  fit  for  the  data 
points  is: 

Y(T)  -  1.3921  -  8.13  x  10"5  T 

This  varies  little  from  a  previous  equation  for  Y  of  propane  and  air  at  ♦  - 
0.9  which  is: 

Y(T)  -  1.3751  -  6.9910  x  10~5  T 

Two  points  were  checked  for  the  unburned  mixture  at  500°K  and  1000°K  for 
isooctane  at  <p  -  1.0.  These  values  of  Y  were  calculated  by  finding  the  slope 
of  h  and  u  at  500°K  and  1000°K  from  using  reference  3  and  using  the 
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relationships: 


These  two  points  are  also  shown  in  figure  2-2. 


2.5  RESIDUAL  GAS 

Residual  gas,  the  exhaust  gas  remaining  in  the  cylinder  after  intake  is 
completed,  affects  the  amount  of  the  new  fuel  and  air  that  can  be  drawn  into 
the  cylinder  and  therefore  affects  the  energy  entering  the  system.  Residual 
gas  fractions  were  calculated  from  the  graphs  of  Toda,  et  al  (4).  Engine  A 
was  chosen  for  its  valve  timing  and  because  it  has  27°  of  valve  overlap  close 
to  the  24°  far  the  Ricardo  engine.  The  compression  ratio  of  8.5  was  chosen 
compared  to  8.41  for  the  Ricardo  engine.  The  residual  gas  fraction  changes 
little  with  the  changing  air  fuel  ratio  so  14.5  was  used  versus  14.37.  The 
final  graph  was  entered  using  the  intake  pressure  and  the  value  of  the 
residual  gas  was  obtained  by  interpolating  for  the  speed.  The  values  for 
residual  gas  are  listed  on  the  engine  data  sheets  (Appendix  A). 

2.6  COMBUSTION  HEAT  RELEASE 

The  net  heat  release  is  obtained  by  considering  the  cylinder  as  a  closed 
system,  modifing  for  the  mass  loss  to  the  crevice,  and  by  considering  the 
combination  of  chemical  energy  added  and  the  heat  transfer  to  the  walls  as  a 
"net  heat  addition".  By  separating  the  heat  transfer  to  the  wall  from  the 
"net  heat  addition",  the  combustion  heat  release  is  obtained: 


* 


The  total  combustion  heat  release  obtained  can  be  compared  with  the  fuel 
energy  added  to  determine  the  validity  of  the  model.  The  best  choice  of 
parameters  can  be  calculated  by  finding  the  error  for  each  of  the  engine  data 
points  and  model  conditions.  The  error  is  found  by  the  following  equation: 

VUIT>'W,  -  <wu  ■ 8 

The  best  fit  is  the  model  condition  for  which  the  square  root  of  the  sum  of 
the  errors  squared  is  the  smallest. 

Best  fit  -  MIN  (Ea2)0,5 

The  model  parameters  varied  are  .  C^,  and  n  in  Woschni's  correlation,  and 
the  crevice  volume  and  the  wall  temperature  in  the  heat  release  analysis. 


CHAPTER  THREE 


ENGINE,  INTRUMEMTATION,  DATA  COLLECTION  AND  PERFORMANCE  VERIFICATION 


3.1  ENGINE  AND  DYNAMOMETER 

The  engine  used  for  data  collection  is  a  Ricardo  HYRDA  MK  ill  (5)  single¬ 
cylinder  spark-ignition  engine  (figures  B-1  and  B-2).  It  is  designed  as  a 
high  speed  research  engine  typical  of  the  engines  used  in  today's  cars.  Its 

O 

compression  ratio  is  8.41 ,  its  displaced  volune  is  496  cor  and  it  has  a  bore 
and  stroke  of  8.57  cm  and  8.6  cm  respectively.  The  engine  has  Bosch  timed 
electronic-fuel-  injection  with  a  separate  electronic  fuel  pimp.  The  fuel 
used  is  Indolene,  a  standard  reference  gasoline.  Indolene  characteristics 
(figure  B-4)  include  a  research  octane  nunber  of  97.2  and  it  is  therefore 
resistant  to  knock.  The  engine  operating  range  is  from  1000  RPM  to  5400  RPM. 
It  has  a  maximun  power  output  of  15  kV. 

The  engine  is  directly  coupled  to  an  EATON  Dynamometer  6000  series 
adjustable-frequency  regenerative  dynamometer  of  100  HP  capacity.  It  has  a 
strain  gage  which  is  calibrated  by  the  dead  weight  method  to  give  a  direct 
readout  of  torque  in  foot  pounds  (lb^-ft).  The  dynamometer  can  operate,. in 
speed  control  mode  which  maintains  a  set  speed  while  allowing  the  load  to 
vary,  or  it  can  operate  in  load  control  while  letting  the  speed  vary.  In  the 
speed  control  mode,  the  dynamometer  can  provide  drive  force  only,  absorb  only, 
or  both  drive  and  absorb  when  In  universal  mode.  Similarly,  the  load  control 
can  provide  either  a  drive  or  an  absorb  load. 


3.2  ENGINE  INSTUMENTATION 

The  pressure  transducer  used  to  measure  cylinder  pressure  is  a  Kistler 
7061  piezoelectric  pressure  transducer  mounted  in  the  cylinder  head  coupled  to 
a  Kistler  charge  amplifier.  The  coupled  set  was  calibrated  by  the  dead  weight 
method  as  per  Lancaster,  et  al  (6).  The  equation  of  the  line  as  calibrated 
for  20  points  from  0  to  1000  psi  is  millivolts  -  1.7  ♦  4.985  x  pressure  (psi) 
with  R2  -  1.0000  (figure  B-5). 

Engine  temperature  measurements  were  made  through  use  of  type  K  chromel 
alumel  thermocouples.  Thermocouples  were  Installed  to  measure  the  oil 
temperature,  the  water  inlet  and  the  water  outlet  temperatures  frcm  the  head 
and  cylinder.  The  air  inlet  and  the  exhaust  temperatures  were  also  measured 
along  with  the  temperature  of  the  fuel  from  the  buret. 

Three  manometers  for  pressure  measurements  were  installed.  These  were: 

(1)  a  mercury  manometer  to  measure  the  time-average  intake  pressure 
Just  downstream  of  the  throttle, 

(2)  a  water  manometer  to  measure  the  air  flow  laminar-flow-element 
pressure  drop  (which  is  upstream  of  the  throttle),  and 

(3)  a  water  manometer  to  measure  the  average  exhaust  pressure. 

The  air  flow  is  measured  by  an  Alcock  viscous  flow  air  meter  which 

measures  the  differential  pressure  across  a  laminar  flow  element.  This 
pressure  difference  is  measured  with  a  special  manometer  that  allows 
variations  in  the  air  flow  ranges.  The  factory  calibration  showed  better  than 

•7 

3 %  linearity.  The  measured  air  flow  in  m  /sec  is  converted  to  gm/sec  at  NTP 
to  obtain  the  measured  air  mass  flow  rate. 
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The  fuel  system  has  a  separate  12  V  fuel  pump  with  a  pressure  regulating 
return  valve  which  maintains  approximately  2  bar  pressure  at  the  fuel 
injection  valve  (figure  B-6) .  The  system  is  set  up  to  draw  and  return  fuel  to 
the  fuel  tank  during  normal  operations  and  while  filling  the  buret.  To 
measure  fuel  flow,  the  supply  and  return  valves  are  simultaneously  switched 
from  the  fuel  tanks  to  the  burets.  The  selection  the  large  or  small  buret  is 
made  by  the  position  of  the  buret  pet cocks.  The  flow  out  of  the  buret  is 
measured  in  milliliters.  The  flow  is  timed  and  is  then  converted  to  a  mass 
flow  rate  by  dividing  by  the  density  of  the  fuel  at  the  measured  temperature 
(figure  B-4  and  Appendix  G). 

A  shaft  encoder  is  connected  to  the  end  of  the  engine  crankshaft  to 
provide  360  pulses  per  revolution  plus  one  larger  pulse  per  revolution.  The 
single  pulse  pa*  revolution  is  set  at  BDC  to  act  as  a  trigger  and  as  a 
reference  point  for  the  data  aquistian  system. 

3.3  DATA  ACQUISTION 

A.  Equipment  hookup.  The  shaft  encoder  and  charge  amplifier  are  easily 
connected  to  the  digitizer.  The  shaft  encoder  output  goes  to  the  tachometer 
interpreter  board  where  the  signal  is  separated  into  two  outputs,  a  pulse  per 
revolution  output  used  as  a  trigger  and  a  crank  angle  degree  output.  The 
pulse  pa*  revolution  output  goes  through  an  amplifier  to  three  places:  the 
ocillosoope  trigger,  to  the  digitizer  stop  trigger  and  to  channel  2  of  the 
digitizer.  The  crank  angle  degree  output  is  connected  to  the  digitizer 
external  clock.  The  pressure  transducer  output  is  connected  to  channel  1  of 
the  ocillosoope  and  to  channel  1  of  the  digitizer. 


The  connection  of  the  cooputor  terminal  and  the  digitizer  to  the  main 
computer  system  is  accomplished  by  the  Sloan  Lab  staff. 

B.  Data  Acquisition  procedures.  After  completing  startup  checks,  the 
engine  is  motored  and  the  throttle  is  set  to  the  desired  load.  From  the 
measured  air  flow,  the  approximate  fuel  setting  is  determined.  Note,  the  air 
flow  drops  about  10  t  between  motoring  and  firing.  The  engine  is  then  fired 
and  the  fuel  and  air  flow  measurements  are  taken.  The  fuel  is  adjusted  until 
the  equivalence  ratio  is  about  1.0.  The  $  (exhaust)  is  calculated  from  the 
measurements  taken  from  the  exhaust  gas  cart  analysis  and  compared  to  $ 
(Intake)  as  measured.  When  these  values  confirm  that  $  «  1.0,  maximum  brake 
timing  (MBT)  is  then  searched  for.  MBT  can  be  found  by  adjusting  the  timing 
to  obtain  the  maximum  measured  brake  torque.  An  alternative  method  for 
determining  MBT  is  from  the  calculated  IMEP  from  /pdV.  After  finding  MBT 
(figure  B-13),  ♦  is  checked  to  ensure  that  it  is  still  equal  to  1.0. 

The  actual  data  acquisition  is  straight  forward.  The  digitizer  parameters 
must  be  verified  and  then  the  digitizer  is  brought  on  line  by  the  ccmputor 
terminal,  the  stop  trigger  is  disconnected  and  data  collection  proceeds.  The 
stop  trigger  is  reconnected  and  the  acquired  data  is  transfered  to  the  VAX 
disc.  The  procedure  is  repeated  for  the  number  of  sets  of  data  required.  The 
system  records  44  cycles  of  data  each  time  and  normally  two  sets  of  44  fired 
cycles  were  collected. 

To  collect  motored  data,  the  fuel  injection  is  turned  off  Just  prior  to 
the  start  of  collection  so  that  the  engine  has  little  time  to  cool  off.  After 
recording  the  motored  data,  the  fuel  injection  is  turned  back  on.  The  $ 
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(Intake)  versus  the  $  (exhaust)  is  cheeked  and  pictures  of  the  pressure  versus 
time  on  the  oscilloscope  are  taken.  Environmental  conditions  and  motored 
conditions  are  checked  and  recorded  and  then  the  engine  is  changed  to  another 
setpoint. 

3.4  ENGINE  SETUP  VERIFICATION  PROCEDURES 

a.  TDC  Location.  The  most  important  item  on  the  engine  requiring 
verification  is  the  top  dead  center  (TDC)  marking.  The  position  of  TDC  is 
checked  by  three  different  methods.  The  simplest  method  is  to  locate  the 
piston  at  the  marked  TDC  and  rotate  the  engine  in  one  degree  increments  about 
TDC  to  +_ 5°.  The  dial  indicator  reading  is  noted  at  each  degree. 

The  second  and  third  methods  are  related.  In  the  second  method,  the 
engine  is  rotated  through  bottom  center  (BC)  to _+  60°  from  TDC.  A  depth 
reading  to  the  piston  crown  is  taken,  then  the  engine  is  rotated  through  BC  to 
61°,  60°,  and  59°  from  TDC.  A  depth  reading  is  taken  at  each  of  these 
positions  and  is  compared  to  the  first  reading  at  60°. 

For  the  third  check,  an  adaptor  was  made  for  the  spark  plug  hole  into 
which  a  sliding  rod  could  be  placed.  The  adaptor  ensured  that  the  rod  is  held 
vertical  and  touches  the  piston  crown  at  the  same  spot  during  the  entire 
dieck.  A  dial  indicator  is  fastened  to  this  rod  and  set  to  about  half  its 
full  scale  reading  with  the  engine  again  positioned  at  _*  60°  from  TDC.  The 
engine  is  rotated  through  BC  until  the  same  reading  is  obtained  on  the  dial 
indicator  on  the  opposite  side  of  TDC.  All  three  checks  were  in  agreement  and 
indicated  that  the  actual  TDC  was  approximately  0.1°  after  the  marked  TDC. 
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This  is  well  within  the  required  accuracy  of  the  data  aquisition  system  and 
the  heat  release  program. 

b.  Phasing.  The  next  item  requiring  verification  is  the  timing  of  the  BC 
pulse  from  the  shaft  encoder.  To  verify  this  timing,  the  pulse  output  is 
connected  to  a  strobetach.  The  engine  is  motored  and  the  BC  pulse  is  checked 
against  the  flywheel  crankangle  markings. 

For  additional  checks  of  the  phasing,  motoring  pressure  data  is  taken. 
From  this  data,  the  average  peak  pressure  location  is  determined.  The  peak 
pressure  location  for  1200  RPM  at  MOT  is  1 .00  +_  0.61  degrees  before  TDC.  The 
log  P-log  V  diagrams  also  were  plotted  from  this  Information  and  the  slopes 
checked  for  curvature  according  to  Lancaster,  et  al  (6).  The  phasing  appears 
correct  from  the  above  checks  (figure  B~7  and  B-8) . 

Another  verification  is  to  check  peak  pressure  versus  PVN  -  constant.  At 
1200  RPM  WOT.  the  transducer  peak  pressure  is  1.2  V  on  the  oscilloscope  which 
equals  240  psig.  240  psig  ♦  15  psi  -  255  psia.  With  Y  -  1.35,  PV*1  yields  a 
value  for  P2  -  260  psia.  These  results  are  very  close. 

c.  and  inlet  pressure.  The  volunetrlc  efficiency  has  also 
checked.  Using  paQ  of  the  air  in  the  engine  cell,  at  2500  WOT  is  80  t 
tar  firing.  The  volumetric  efficiency  for  motoring  is  typically  10  t  higher 
due  to  10  %  higher  air  flow  during  motoring. 

Prior  to  taking  data,  the  data  acquisition  system  requires  a  correction 
for  bias  and  scaling.  The  scaling  converts  the  charge  amplifier  input  though 
the  digitizer  into  a  value  of  BITS/VOLT  for  use  in  the  data  acquisition 
system.  The  system  is  scaled  as  follows: 


4.985  m  V  _  14.69595  psi  ..  4096  bits  ^  _  bits 

psl  x  atm  x  10  Volts  *  JU‘0Uf  atm" 

The  data  acquistion  requires  a  negative  scaling  factor  in  atm/bits.  So  the 
scaling  factor  is: 


-  3.33255  x  10-2 


atm 

bits 


The  bias  is  set  so  that  4096  bits  equals  zero  atmospheric  pressure  (zero 
Volts  output  from  the  charge  amplifier).  The  full  scale  range  of  the  data 
acquisition  system  is: 

4096  bits  x  3.33255  x  IQ-2  -  136.5014  atm 


Therefore  the  bias  is  set  at  136.5  atm.  Figure  8-8  shows  a  log  P-log  V 
diagram  drawn  prior  to  setting  the  scale  and  bias.  The  curved  compression  and 
expansion  lines  indicate  that  an  incorrect  intake  pressure  was  used  for  this 
graph.  Figure  shows  the  same  data  plotted  after  correcting  for  scale  and 
bias. 

The  intake  pressure  measured  on  the  intake  manometer  is  not  the  actual 
pressure  during  the  engine  intake  process.  It  is  the  average  pressure  in  the 
Intake  port  which  equals  Pj  from  about  12°  before  TDC  to  56°  after  BC  (i.e. 
248°  of  720°)  and  approximately  equals  atmospheric  pressure  the  rest  of  the 
cycle.  The  initial  approximation  for  P^  is  made  using  the  assumption  that  the 
air  flow  is  proportional  to  the  pressure.  Pj  is  approximately  as  follows: 


Px (During  intake)  -  PI(WT) 


a  (WOT) 
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Using  this  value  for  P^,  the  log  P-log  V  diagrams  for  the  fired  data  are 

graphed  with  Pj.  varied  by  5  %  increments  of  its  estimated  value.  The  slopes 

of  the  compression  and  expansion  slope  are  checked  for  curvature  and  the 

intake  and  expansion  lines  are  checked  to  ensure  that  the  pressure  in  the 

cylinder  is  less  than  during  intake  and  greater  than  patm  during  exhaust 

as  per  Lancaster,  et  al  (6).  The  values  that  provide  the  best  fit  log  P-log  V 

diagrams  are  the  values  used  for  P  •  in  the  heat  release  analysis. 

d.  Torque  and  MEP.  Verification  of  the  engine  performance  can  be  made  by 

comparing  the  torque  readout  from  the  dynamometer  with  the  engine  IMEP. 

Defining  the  terms  IMEP _ as  the  IMEP  calculated  from  the  pressure  data  for 

gross 

the  compression  and  expansion  strokes  and  IMEPnet  as  the  calculated  IMEP  for 
the  entire  engine  cycle,  the  difference,  IMEP  -  IMEP ,  equals  the 

nCU 

pumping  MEP  or  PMEP.  The  IMEP _ calculated  from  the  /pdV  can  be  compared 

gross 

to  the  IMEP _  calculated  from  the  sun  of  the  brake  torque,  the  friction, 

gross 

and  the  punping  MEP.  These  IMEP^^  values,  when  plotted  against  each  other 
should  result  in  a  straight  line  with  a  slope  of  t  (figure  B-9).  The 
IMEPgross  from  the  is  calculated  by  integrating  the  measured  pressure 
over  the  known  volume  change.  The  IMEP  (from  torque)  is  calculated  by  the 
following  method: 

Brake  MEP  (from  motored  torque)  ■  Rubbing  friction  MEP  +  PMEP  (motored) 
Therefore: 

Rubbing  friction  MEP  -  Brake  MEP  (motored  torque)  -  PMEP  (motored) 
IMEPgroag  (ft*®*  torque)  -  Rubbing  friction  MEP  ♦  Brake  MEP  (from  fired 
torque)  +  PMEP  (fired) 
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The  result  of  the  calculation  is  that  the  IMEP  (torque)  is  less  than  the 
IMEP(/pdV)  at  lower  values  of  IMEP,  but  are  approximately  equal  at  higher 
values  of  IMEP. 

IMEP  calculated  from  pdV  and  the  BMEP  calculated  from  the  torque  can 
gross 

be  plotted  against  the  air  flow  normalized  by  the  WOT  air  flow  (ma/*al_T) . 

The  linear  least  squares  curve  fit  for  IMEP _ is: 

gross 


A 

-4.71  ♦  152.8  [  ] 

“a  WOT 


and  for  the  EMEP  from  the  measured  torque  is: 

m 

BMEP  (torque)  -  -27.19  ♦  144.79  [  — - ] 

*a  WOT 

The  slopes  of  the  plots  of  IMEP _ _  and  BMEP  (torque)  from  the  least  squares 

gross 

fit  are  nearly  equal  (152.8  versus  144.79).  These  results  indicate  that  the 
calculated  IMEP  from  the  pressure  data  and  the  IMEP  from  the  torque  are 
compatible  and  that  the  motoring  IMEP  (friction  and  pumping)  is  approximately 
constant  over  the  engine  load  and  speed  set  points  only  increasing  slightly 
(from  25  to  30  psi)  with  load  (figure  B-10). 

e.  F/A  measurements  and  $.  An  accurate  measurement  of  the  fuel  flow  and  of 
the  air  flow  is  required  during  data  acquisition  to  ensure  reasonable  results 
are  obtained  for  the  heat  release  analysis.  The  first  step  is  to  determine 
fuel  flow  versus  fuel  dial  settings  at  different  speeds.  Initial  plots  were 
made  by  running  the  engine,  injecting  the  fuel  into  a  clean  container,  and 


measuring  the  fuel  flow  rate  at  various  speed  and  fuel  flow  settings.  Then  as 
data  is  collected,  additional  fuel  flow  rates  at  different  settings  can  be 
collected.  The  flow  rate  chart  (figure  G-3, Appendix  G)  can  be  used  as  an 
approximate  value  at  which  to  set  the  fuel  for  a  desired  flow  rate. 

The  fuel  flow  rate  is  the  basis  for  heat  addition  into  the  engine: 
heat  addition  -  (LHV) .  It  is  therefore  important  to  ensure  that  the 
measurement  of  the  fuel  flow  rates  is  as  accurate  as  possible.  Part  of  the 
calculation  for  mass  flow  rate  is  to  divide  the  volune  flow  rate  by  the 
density  of  the  fuel  at  the  fuel  temperature.  The  density  of  the  fuel  is 
calculated  by  using  linear  curve  for  the  change  of  specific  gravity  as  a 
function  of  temperature  as  suggested  by  Bolt,  et  al  (7)  and  applying  that 
curve  to  the  measured  API  of  the  fuel. 

The  resulting  equation  is  density  of  lndolene: 

density  •  0.9990[0.7523  ♦  (15.6  -  fuel  temperature(°C))0. 00081] 

A  measured  quantity  of  fuel  was  weighed  to  check  the  accuracy  of  the  density 
calculation.  At  20°C,  the  calculated  density  of  indolene  is  .7480  g/cm^  and 

•3 

the  measured  density  is:  .7455  g/cm.  The  error  is  les3  than  0.4  % 

The  air  flow  rate  is  calculated  from  a  measurement  by  a  manometer.  The 

m 

fuel  flow  rate  is  measured  from  a  buret  graduated  in  millileters.  The  ambiant 
air  pressure  affects  the  first  while  the  temperature  affects  both 
measurements.  The  equivalence  ratio  can  be  calculated  from  the  measured 
flews. 

Verification  of  the  measured  flow  ratio  is  completed  by  calculating  the 
equivalence  ratio  based  on  the  exhaust  gas  anaylsis.  After  the  engine  has 
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stabilized  at  a  set  point,  the  fuel  and  air  flow  rates  should  be  measured 
while  simultaneously  measuring  the  exhaust  gas  constituents.  Fran  these 
measured  values,  the  equivalence  ratio  is  calculated  and  compared.  The  error 
should  be  less  than  +_  3  $  •  The  equivalence  ratio  calculated  from  the  measured 
input  flow  versus  the  exhaust  gas  cart  anaylsis  equivalence  ratio  have  been 
plotted  (figure  B— 1 1 ) -  The  calculated  exhaust  $  compares  favorably  with  the 
calculated  input  $  for  $  less  than  1.1. 

The  accuracy  of  the  exhaust  gas  analysis  equipment  is  checked  by  plotting 
the  individual  exhaust  species  concentration  versus  the  ideal  species 
concentration  as  a  function  of  the  exhaust  equivalence  ratio  (figure  B-12). 

The  equation  used  to  calculate  the  ideal  emissions  at  the  equivalence  ratio  <t> 
is: 

V,1.83  *  *£75  CO,.  3.76  «2)  _> 

♦ 

a  C02  +  b  H20  +  c  C02  +  d  H2  +  37. W  \ 

♦ 


The  ratio  of  carbon  monoxide  to  hydrogen  was  assumed  to  be  three:  C  -  3d. 
Carbon  and  hydrogen  balance  results  in  these  equations: 
a  +  c  -  7 


2b  ♦  2d  -  11.83 


9.9575 

♦ 


2a  +  b  ♦  c 


Solving  these  last  four  equations  simultaneously  at  a  chosen  <p  yields  the 
ideal  emissions  at  that  equivalence  ratio. 

A  more  complete  check  is  to  use  the  ratio: 
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[co]  [h2o] 

[h2T  [coTT  “3‘5 

to  compute  the  ideal  emission  versus  Also  a  calculation  can  be  made  where 
the  measured  carbon  monoxide  and  the  measured  oxygen  can  be  combined  and 
compared  to  the  ideal  emissions  by  a  carbon  atom  count.  These  last  two  checks 
were  not  completed  for  this  analysis. 
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CHAPTER  FOUR 

ANALYSIS,  RESULTS,  AND  DISCUSSION 


4.1  COMPUTER  PROGRAM 

The  computer  program  used  for  the  heat  release  analysis  is  the  computer 
program  called  QCAL  at  the  MIT  Sloan  Automotive  Lab  started  by  A.Giovanettl(8) 
and  developed  by  J.  Gatowski.  The  program  was  developed  using  a  one  zone 
model  for  the  gas  in  the  cylinder.  It  has  been  developed  to  be  compatible 
with  the  SLN  data  acquisition  program  at  MIT  Sloan  Automotive  Lab. 

The  heat  release  program  calculates  the  mass  of  gas  introduced  into  the 
cylinder.  It  smooths  the  pressure  data  curve  and  calculates  the  differential 
pressure  change  (dp)  for  each  time  step.  It  provides  the  choice  of  individual 
cycle  analysis  or  of  average  cycle  analysis.  It  calculates  the  instantaneous 
heat  transfer  to  the  walls,  the  work  done  in  the  cylinder  against  the  piston 
and  corrects  for  the  flow  into  the  crevice  and  for  the  crevice  heat  loss.  It 
then  integrates  these  values  to  find  the  combustion  heat  release  as  detailed 
in  sections  2-1  and  2-6. 

The  computer  program  requires  the  motored  engine  data  file  and  the  fired 
engine  data  file  as  Inputs.  It  prints  a  stannary  file  with  a  list  of  engine 
data,  environmental  conditions  and  heat  release  inputs  that  can  be  modified 
for  the  present  run.  The  output  of  the  program  is: 

(1) a  printed  stannary  listing  the  inputs  used  and  the  heat  release  values 
calculated, 

(2)  a  plot  file  of  the  combustion  heat  release,  and 


(3)  a  plot  file  of  the  energy  in  the  cylinder  modified  by  the  mass 
effects. 

4.2  CALIBRATION  OF  THE  MODEL 

1 .Variations  of  the  Model  Parameters.  There  are  seven  variables  that  are 

required  as  inputs  to  the  heat  release  program.  These  are:  C1 ,  C2»  V^,  Vc, 

Y,  Tw,  and  the  exponent  (n)  in  the  heat  transfer  equation  (section  2.3).  Two 

of  these,  V  and  Y,  can  be  calculated  prior  to  the  heat  release  analysis.  V 
c  c 

can  be  calculated  from  the  engine  dimensions  or  can  be  approximated  from  the 
compression  ratio  by  using  the  relationships  VQ  ♦  VQ  -  VT  and  VT/Vc  -  rQ  to 
get: 


where  V  is  the  clearance  volume,  V_  the  displacement,  V_  the  volume  at  BC  and 
c  u  1 

r  the  compression  ratio  (equal  to  8.41  :  1).  Y  is  calculated  as  shewn  in 
0 

section  2.4. 

a.  Changing  C^  Increasing  C^  causes  the  heat  transfer  to  increase  as  Q  * 
C^hAAT.  If  the  heat  transfer  is  estimated  to  be  larger  than  the  actual  heat 
transfer,  then  the  combustion  heat  release  will  continue  to  rise  after 
combustion  has  been  completed.  Therefore,  the  end  slope  of  the  combustion 
heat  release  profile  increases  (figure  C-1  Appendix  C). 

b.  Changing  C^  Increasing  C2  increases  the  heat  transfer  during 
combustion  and  expansion.  A  greater  heat  transfer  to  the  walls  during 
combustion  forces  the  combustion  heat  release  during  combustion  to  increase. 
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The  result  is  a  large  increase  in  the  total  value  of  combustion  heat  release 
(figure  C-2  Appendix  C). 

c.  Changing  Vcr  The  effect  of  crevice  volune  as  detailed  in  section  2.2 
is  both  mass  and  energy  loss  from  the  control  volune.  The  increased  loss  to 
the  crevice  due  to  an  increase  in  the  crevice  size,  results  in  a  small 
Increase  of  the  value  of  the  combustion  heat  release  over  the  entire 
combustion  heat  release  curve  and  a  drop  in  the  total  energy  in  the  cylinder, 
especially  during  the  high  pressure  portion  of  the  engine  cycle  (figure  03 
Appendix  C). 

d.  Changing  T^  Changing  the  wall  temperature  over  a  plausible  range  has  a 
small  effect  on  the  combustion  heat  release.  As  Ty  increases,  the  temperature 
difference  between  the  gas  and  the  walls  decrease  resulting  in  a  decrease  in 
the  heat  transfer  after  combustion.  An  increase  in  T^  causes  a  small  decrease 
in  the  slope  at  the  end  of  the  combustion  heat  release  curve  (figure  C-*» 
Appendix  C). 

e.  Changing  the  exponent  Woschni's  heat  transfer  correlation  was 
rewritten  in  section  2.3  to  incorporate  an  exponent  (n)  that  could  be  varied 
while  maintaining  the  relationships  between  B,  P,  T  and  W.  The  affect  of 
changing  n  is  opposite  to  the  affect  of  changing  both  C^  and  C2>  Decreasing  n 
increases  the  heat  transfer,  thus  the  combustion  heat  release  curve  displays 
an  increase  in  heat  release  during  combustion  and  an  increase  after  combustion 
(figure  C-5  Appendix  C). 

The  reference  pressure  used  is  the  intake  pressure.  The  intake  pressure 
has  only  a  small  affect  on  the  combustion  heat  release  curve  when  the  pressure 
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used  in  calculations  is  near  the  correct  intake  pressure.  However,  a  large 
error  in  the  intake  pressure  presents  itself  as  a  change  in  the  slope  of  the 
combustion  heat  release  curve  before  and  after  combustion.  As  the  value  for 
PT  used  in  the  analysis  increases,  the  slope  before  combustion  decreases  and 
the  slope  after  combustion  increases  (figure  C-6  Appendix  C). 

2.  Best  Fit  of  the  Model  Parameters  The  data  for  one  engine  operating 
set  point  was  initially  chosen  for  testing  in  the  heat  release  analysis  to 
determine  the  correct  model  parameters.  The  original  setpoint  chosen  was  1500 
RPM  with  Pj-0.35  atm.  The  parameters  listed  in  section  4.2.1  were  varied  in 
the  order  detailed  there  to  determine  the  set  of  parameters  that  yielded  the 
desired  combustion  heat  release  profile.  The  profile  sought  was  flat  with 
zero  value  before  combustion,  nearly  flat  after  combustion,  and  approached  but 
did  not  exceed  the  energy  in  the  cylinder  curve.  When  a  set  of  model 
parameters  was  found  that  sucessfully  modeled  the  data  at  1500  RPM  with 
Pj-0.35  atm,  the  chosen  parameters  were  then  checked  against  the  other  five 
sets  of  data.  The  model  parameters  were  individually  varied  for  each  set  of 
engine  data  to  verify  the  effects  of  varying  model  paramters  for  that  engine 
data  set.  The  results  of  varying  the  parameters  were  found  to  be  similar  for 
each  set  of  engine  data. 

The  original  set  of  model  paramters  sucessfully  modeled  the  heat  release 
for  all  six  sets  of  engine  data  as  determined  by  a  visual  qualitative  check  of 
the  graphs.  To  check  if  another  set  of  model  parameters  would  provide  a 
better  qualitative  fit,  the  set  of  chosen  parameters  were  varied  and  the  six 
sets  of  data  again  analyized  with  the  different  variables.  The  data  was 


compared  by  speed  and  load  as  displayed  in  the  graphs  in  Appendix  F.  The  new 
set  of  graphs  obtained  from  changing  the  model  parameters  was  compared  to  the 
original  set  of  graphs  to  determine  which  model  provided  a  better  fit.  The 
first  chosen  sets  of  parameters  continued  to  provide  the  best  qualitative  fit. 

A  quantitative  fit  was  desired  (Appendix  E).  The  fit  chosen  was  to 
compare  the  maximum  value  of  the  combustion  heat  release  (determined  from  the 
heat  release  analysis)  with  the  chemical  energy  released  value  (mass  of  fuel 
multiplied  by  the  lower  heating  value  multiplied  by  the  combustion 
efficiency).  The  difference  between  the  two  values  equals  the  error,  e,  which 
was  normalized  by  the  chemical  energy  released  value  for  each  curve  to  remove 
the  effects  of  speed  and  load.  A  least  squares  error  method  was  used  to 
determine  the  error  for  each  set  of  model  parameters  for  the  six  data  sets. 
This  was  done  by  summing  the  values  of  the  error  squared  for  each  of  engine 
data  sets  and  then  taking  the  square  root  of  this  summation.  The  set  of  model 
parameters  which  yields  the  smallest  error  for  the  range  of  the  engine  data, 
is  the  set  of  parameters  that  give  the  "best  fit".  Using  the  best  fit  set  of 
model  parameters  results  in  the  combustion  heat  release  curves  which  provide 
the  closest  approach  to  the  chemical  energy  released  lines. 

4.3  RESULTS 

Graphs  displaying  the  results  of  the  heat  release  analysis  for  each  of  the 
engine  operating  points  are  presented  in  Appendix  D.  There  are  three  lines  on 
each  of  these  graphs. 

The  first  of  these,  which  starts  near  the  top  on  the  left  edge  and  runs 


horizontal  with  a  dip  in  the  middle,  is  a  line  representing  the  energy  in  the 
cylinder,  the  "cylinder  energy"  line.  This  line  is  the  value  of  the 
instantaneous  mass  of  fuel  in  the  cylinder  multiplied  by  the  lower  heating 
value  of  the  fuel.  This  line  starts  out  equal  to  the  maximum  fuel  energy 
introduced  into  the  cylinder,  but  then  dips  as  the  pressure  rises  inside  the 
cylinder  and  forces  mass  (and  therefore  energy)  out  of  the  cylinder  control 
volume  and  into  the  crevice  (section  2.2).  The  value  for  the  line  rises  in 
the  later  part  of  the  curve  as  mass  returns  to  the  cylinder  volune  from  the 
crevices.  The  cylinder  energy  line  in  the  middle  portion  of  the  graph  where 
it  dips,  represents  the  mlnlnun  energy  inside  the  cylinder.  This  is  because 
some  of  the  mass  forced  into  the  crevice  may  be  burned  gas  and  therefore 
removes  less  energy  than  unburned  mass  when  it  leaves  the  cylinder.  The  right 
hand  end  point  should  not  be  greater  than  the  left  hand  end  point  but  may  not 
return  to  the  initial  left  hand  aid  value  due  to  some  mass  remaining  in  the 
crevice  during  expansion  and  due  to  the  crevice  heat  loss  affects  (section 
2.2). 

The  second  line  on  the  graph,  the  short  horizontal  line  near  the  top  right 
hand  side,  is  the  "chemical  energy  released"  line.  This  line  represents  the 
maximum  value  that  the  combustion  heat  release  curve  is  expected  to  achieve. 
The  value  for  this  line  is  the  value  of  the  fuel  energy  introduced  into  the 
cylinder  (mass  of  fuel  multiplied  by  the  lower  heating  value  of  fuel) 
multiplied  by  the  combustion  efficiency  as  calculated  from  the  exhaust  gas 
analysis.  The  combustion  efficiency  is  computed  by  calculating  the  chemical 
energy  remaining  in  the  exhaust  components  of  hydrogen,  carbon  monoxide,  and 


hydrocarbons,  and  from  these  values  determine  the  fraction  of  the  chemical 
energy  introduced  into  the  cylinder  that  was  burned.  Values  for  the 
combustion  efficiency,  the  fuel  energy  introduced  into  the  cylinder,  and  the 
chemical  energy  released  line  are  listed  in  the  engine  data  table  (Appendix 
A).  The  line  for  the  chemical  energy  released  is  the  line  that  the  combustion 
heat  release  curve  should  asymptotically  approach. 

The  third  line  is  the  combustion  heat  release  profile.  This  is  the  line 
representing  the  accumulative  quantity  of  energy  released  by  combustion.  This 
line  is  expected  to  be  flat  and  equivalent  to  zero  until  combustion  starts. 

It  is  then  expected  to  have  a  rapid  rise  over  30-40  degrees  of  crank  angle 
after  which  it  should  again  flatten  out  with  a  slope  of  zero  and  a  value  equal 
to  the  chemical  energy  released  line.  Oxidation  can  occur  after  combustion 
and  into  the  exhaust  process  so  that  the  combustion  heat  released  curve  can 
have  a  positive  end  slope.  Recombination  of  exhaust  products  can  occur  which 
would  allow  the  combustion  heat  release  curve  to  go  above  the  chemical  energy 
released  line  and  allow  it  (the  combustion  heat  released  curve)  to  have  a 
negative  aid  slope.  It  should  then  asymptotically  approach  the  chemical 
energy  released  line  from  above.  The  combustion  heat  release  curve  cannot  be 
greater  than  the  total  fuel  energy  introduced  into  the  cylinder.  Ideally,  it 
should  not  go  above  either  the  cylinder  energy  curve  or  the  chemical  energy 
released  line  and  should  asymptotically  approach  the  chemical  energy  released 
line  as  the  combustion  process  completes. 

a.  1000  RPM  Pj-0.36  atm  (Figure  D-1B  Appendix  D).  The  combustion  heat 
release  curve  reaches  96  %  of  the  value  of  the  chemical  energy  released  and 
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9 2$  of  the  cylinder  energy  line.  The  slope  after  combustion  is  very  flat  with 
a  very  small  turndown  of  the  slope  at  the  end  of  the  combustion  heat  release 
curve.  This  turndown  diverges  from  the  chemical  energy  released  line. 

b.  1500  RPM  P^-0.35  atm  (Figure  D-2B  Appendix  D).  The  combustion  heat 
release  curve  reaches  103  %  of  the  fuel  energy  released  curve,  94  %  of  the 
cylinder  energy  curve,  and  has  a  positive  slope  after  combustion  is  complete. 
This  positve  slope  diverges  from  the  chemical  energy  released  line  but  the 
combustion  heat  release  remains  within  4  %  at  the  fuel  energy  released  line. 

c.  1500  RPM  Pj-,44  atm  (Figure  D-38  Appendix  D).  The  combustion  heat 
release  reaches  99  %  of  the  chemical  energy  released  line,  92  %  of  the 
cylinder  energy  curve,  and  it  has  a  positive'  slope  so  that  it  continues  to 
approach  the  chemical  energy  released  line.  A  dip  appears  in  the  combustion 
heat  release  curve  just  after  combustion  is  completed. 

d.  1500  RPM  WOT  (Figure  D-4B  Appendix  D).  The  combustion  heat  release 
curve  for  WOT  displays  the  same  general  trends  as  does  the  1500  RPM  curve  with 
Pj-,44  atm.  It  reaches  98  %  of  the  chemical  energy  curve,  94  %  of  the 
cylinder  energy  curve,  and  has  the  same  dip  Just  after  completion  of 
combustion.  During  the  final  portion  of  the  combustion  heat  release  curve  the 
slope  becomes  negative,  diverging  from  but  remaining  within  4  %  of  the 
chemical  energy  released. 

e.  2500  RPM  P^^.35  atm  (Figure  D-5B  Appendix  D).  The  combustion  heat 
release  curve  reaches  104  %  of  the  chemical  energy  released  and  95  %  of  the 
cylinder  energy.  The  overall  shape  of  the  curve  is  well  modeled.  The  slope  of 
the  combustion  heat  release  curve  levels  off  sifter  combustion  and  tends  to 


climb  slightly  at  the  end.  The  slope  before  combustion  is  equal  to  zero  as  is 
expected  for  the  model,  but  is  offset  by  approximately  5  Joules.  This 
phenomena  was  only  displayed  for  this  set  of  engine  data. 

f .  2500  RPM  WOT  (Figure  D-6B  Appendix  D) .  The  combustion  heat  release 
curve  reaches  102%  of  the  chemical  energy  released  line  and  contacts  the  curve 
for  energy  in  the  cylinder  before  it  levels  off  at  98  %  of  the  cylinder 
energy.  Near  the  end  of  the  curve,  the  slope  is  negative  and  the  combustion 
heat  release  curve  approaches  the  chemical  energy  released  line.  The  curve  is 
well  modeled  except  that  the  slope  is  greater  than  zero  prior  to  combustion. 
Two  possiblities  indicated  are  that  either  too  small  an  inlet  pressure  was 
used  or  too  large  a  crevice  volume  was  used. 

4.4  COMPARISONS 

The  individual  curves  can  be  analyized  in  different  ways.  The  curves  of 
the  energy  in  the  cylinder,  the  chemical  energy  release,  and  the  combustion 
heat  release  are  compared  versus  (Slanging  speed  at  two  loads:  (a)  with  inlet 
pressure  at  approximately  0.35  atmosphere  and  (b)  at  WOT.  The  curves  can  also 
be  compared  at  steady  speed  with  varying  looads  at  1500  RPM  and  2500  RPM. 

The  comparisons  are  made  by  observing  the  similiarities  in  the  slope  of 
the  curves  of  the  combustion  heat  release  as  either  the  speed  or  the  load 
changes.  The  shapes  should  be  the  same  with  similar  rates  of  combustion  heat 
release  both  during  combustion  and  after  combustion.  The  maximun  value  of 
each  heat  release  curve  should  reach  approximately  the  same  percentage  of  its 
chemical  energy  released  line.  The  curves  can  also  be  compared  relative  to 


the  percentage  of  the  cylinder  energy  that  each  achieves.  In  general,  each 
individual  curve  is  well  modeled  by  the  heat  release  analysis  and  the  curves 
scale  well  with  load  and  speed. 

a.  Varying  speed  at  P^^O.35  atm  (Figure  F-1  Appendix  F)  As  can  be  seen  in 
the  figure,  the  shape  of  the  curves  is  very  similar.  All  three  combustion 
heat  release  curves  show  a  similar  rate  in  the  rise  of  the  combustion  heat 
release  and  have  final  slopes  that  are  nearly  equivalent.  The  combustion  heat 
release  curve  at  1500  RPM  and  2500  RPM  match  even  better.  They  are  nearly 
identical  and  both  achieve  104  %  of  the  chemical  energy  released.  They  also 
both  reach  94  %  of  the  cylinder  energy.  The  scaling  for  speed  between  1500 
RPM  and  2500  RPM  is  ideal  in  this  case. 

It  is  obvious  that  the  timing  of  the  1000  RPM  line  was  10-15  degrees  of 
crank  angle  later  than  opt lmm.  The  shape  of  the  combustion  heat  release  and 
the  cylinder  energy  curves  are  similar  to  those  for  1500  RPM  and  2500  RPM  but 
have  10-15  degrees  phase  shift.  The  1000  RPM  heat  release  curve  achieved  a 
smaller  percentage  of  the  chemical  energy  released  and  a  smaller  percentage  of 
the  cylinder  energy  value  than  did  the  1500  RPM  and  2500  RPM  curves.  The 
difference  for  1000  RPM  is  that  <p  equals  0.96  with  n  (combustion)  equal  to  96 
%  compared  to  1.03  and  91  %  respectively  for  the  1500  RPM  and  2500  RPM  curves. 

b.  Varying  speed  at  WOT  (Figure  F-2  Appendix  F)  The  rate  of  combustion 
heat  release  is  identical  for  both  1500  RPM  and  2500  RPM.  Both  have  a 
positive  slope  after  combustion  which  changes  to  a  negative  slope  after 
approximately  100°  of  crank  angle.  The  energy  in  the  cylinder  and  the  chemical 
energy  released  lines  increase  with  speed,  but  the  increase  in  the  combustion 
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heat  release  is  greater  than  the  increase  in  the  energy  supplied  (a  7  % 
increase  of  combustion  heat  release  versus  a  3  >  Increase  in  fuel  energy). 

The  combustion  efficiencies  for  both  cases  are  identical  but  the  equivalence 
ratio  for  1500  RPM  is  2  $  greater. 

While  the  overall  scaling  for  speed  is  generally  good,  there  is  a  4  % 
difference  at  WOT.  One  parameter  that  could  be  modified  to  reduce  this 
difference  is  the  wall  temperature,  which  is  expected  to  increase  with  speed. 
Modifying  the  wall  temperature  for  each  speed  would  lower  the  combustion  heat 
release  at  2500  RPM  and  improve  the  speed  scaling. 

c.  Varying  load  at  1500  RPM  (Figure  F-3  Appendix  F)  The  combustion  heat 
release  curves  model  well  for  changing  load.  All  these  curves  start  their 
initial  rise  at  approximately  the  same  crank  angle.  The  higher  loads  rise 
faster  to  achieve  the  greater  combustion  heat  release  in  approximately  the 
same  number  of  crank  angle  degrees  as  the  lighter  loads.  They  all  have 
positive  slopes  after  combustion.  The  dip  in  the  combustion  heat  release 
displayed  inmediately  after  combustion  completes  on  the  WOT  curve,  can  be  seen 
to  a  lesser  extent  at  Pj-.44  atm.  All  three  curves  achieve  approximately  93  % 
of  the  cylinder  energy.  Due  to  the  increasing  combustion  efficiency  with  the 
increasing  load,  the  values  of  the  chemical  energy  released  rise  faster  than 
the  combustion  heat  release.  The  value  for  each  of  the  combustion  heat 
release  varies  less  than  4  %  from  its  respective  chemical  energy  released 
line. 

d.  Varying  load  at  2500  RPM  (Figure  F-4  Appendix  F)  The  two  combustion 
heat  release  curves  at  Pj-.35  and  at  WOT,  also  look  similar  in  appearance. 


They  both  have  a  alight  rise  after  combustion  and  aid  with  a  small  negative 
slope.  Unlike  the  curves  at  1500  RPM,  both  combustion  heat  release  curves  at 
2500  RPM  extend  above  the  chemical  energy  released  line,  but  both  decrease 
near  the  end  to  reapproach  the  chemical  energy  released  line.  The  final 
values  for  the  combustion  heat  release  curves  are  within  2  %  of  the  chemical 
energy  released  line. 

e.  General  comments  (Figure  F-5  Appendix  F)  The  heat  release  analysis 
successfully  models  the  combustion  heat  release  for  the  range  of  engine  data 
checked.  The  combustion  heat  release  curves  are  all  flat  before  combustion, 
then  display  a  large  rise  for  combustion,  and  then  end  with  a  slope 
approximately  equal  to  zero.  The  combustion  heat  release  curves  reach  92  to 
94  {  of  the  energy  introduced  into  the  cylinder  and  are  within  4  %  of  the 
chemical  energy  released  line. 

The  scaling  with  load  is  very  good.  The  shape  of  the  curves  at  each  load 
is  mirrored  at  the  other  loads  for  that  speed.  The  percentage  of  the  energy 
Introduced  into  the  cylinder  that  appears  in  the  combustion  heat  release  is 
consistant  throughout  the  testing  load  range. 

The  scaling  with  speeed  is  very  good  at  low  load  with  each  curve  reflecting 
the  other  curves  at  that  load.  At  WOT,  the  increase  in  combustion  heat 
release  is  greater  than  the  increase  of  the  energy  introduced  into  the 
cylinder.  Varying  the  wall  temperature  or  varying  the  exponent  (in  equation  1 
section  2.3)  may  fine  tune  the  scaling  with  speed. 
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CHAPTER  FIVE 
SUMMARY  AND  CONCLUSION 


Engine  performance  and  the  details  of  the  engine  combustion  process  can  be 
studied  by  using  measured  cylinder  pressure  at  known  volumes  in  a  heat  release 
analysis.  The  heat  release  analysis  uses  the  First  Law  of  Thermodynamics  and 
the  ideal  gas  law  to  determine  the  apparent  heat  release  in  the  cylinder. 
Modifying  the  apparent  heat  release  for  heat  transfer  to  the  walls  and  for  the 
effect  of  mass  and  energy  flows  into  the  crevices,  the  heat  release  from 
combustion  can  be  determined.  The  equation  for  the  combustion  heat  release 
from  section  2.6  is: 

V* 

cr 

where  6^  is  the  modification  for  the  heat  transfer  to  the  walls  using 
Woshni's  (1)  equation  the  combination  of  the  second  and  third  terms  is  the 
work  done  in  the  cylinder  against  the  piston,  and  the  last  term  is  the 
modification  for  the  crevice  effect. 

The  Ricardo  MK  III  engine  was  instrumented  for  data  collection  and 
prepared  for  operation.  The  engine  was  operated  and  the  operating  parameters 
were  (Sleeked  for  consistancy  with  satisfactory  results.  Engine  cylinder 
pressure  was  measured  for  both  motored  and  fired  conditions  at  six  operating 
points  in  the  range  of  1000  RPM  to  2500  RPM  and  in  the  load  range  from  intake 
pressure  equal  to  0.35  atmosphere  to  WOT.  Forty  four  cylces  of  engine 
pressure  data  were  averaged  for  each  engine  operating  point  and  the  average 
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pressure  data  was  used  in  the  Sloan  Automotive  Laboratory  heat  release  model 
to  determine  the  combustion  heat  release  profiles.  The  model  parameters 
listed  in  section  4.2.1  were  scaled  to  fit  the  engine  data.  The  value  C^, 
used  to  modify  the  heat  transfer  relation,  was  found  to  be  of  the  order  of 
unity  as  expected. 

The  combustion  heat  release  was  successfully  modeled  with  one  set  of  model 
parameters  for  the  range  of  the  data  checked.  The  resulting  profiles  were 
flat  with  a  zero  value  before  combustion,  they  had  a  rapid  rise  in  heat 
release  during  combustion,  and  they  leveled  off  after  combustion  at 
approximately  93  t  of  the  value  of  the  fuel  energy  introduced  into  the 
cylinder.  The  criteria  used  to  select  the  best  set  of  model  parameters  was  a 
least  squares  error  fit  to  the  value  of  the  chemical  energy  released.  The 
chemical  energy  released  was  determined  from  the  combustion  inefficiency 
(calculated  from  the  chemical  energy  of  the  measured  hydrogen,  carbon  monoxide 
and  hydrocarbons  remaining  in  the  exhaust  gas)  by  multiplying  the  fuel 
introduced  into  the  cylinder  by  the  lower  heating  value  of  the  fuel  and  by  the 
combustion  efficiency  (1 .0  -  inefficiency) .  The  ideal  combustion  heat  release 
should  asymptotically  approach  the  value  of  the  chemical  energy  released.  The 
selected  model  parameters  provided  heat  release  profiles  that  were  within  4  % 
of  the  chemical  energy  released. 

The  heat  release  model  is  simple  to  use  and  provides  accurate  results  for 
this  engine.  To  use  this  model  for  another  engine,  the  model  parameters 


should  be  rescaled 
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APPENDIX  A 
ENGINE  DATA  SUfrMARY 


Speed  (RPM) 

1000 

1500 

1500 

1500 

2500 

2500 

Intake  pressure  (atm) 

0.36 

0.35 

0.44 

WOT 

0.35 

WOT 

Spark  Timing  0  BTDC 

30 

45 

40 

*22 

65 

30 

Torque  lb  -ft 
brake/motored 

3.2 

-TPT 

— 

6.0 

-5.2 

23.6 

^C4 

3.2 

-6.6 

— 

IMEP  gross  (atm) 

2.67 

2.86 

3.67 

9.75 

2.91 

10.42 

Fuel  Flow  (g/s) 

0.082 

0.126 

0.163 

0.414 

0.215 

0.710 

Air  Flow  (g/s) 

1.131 

1.724 

2.300 

5.787 

2.968 

9.927 

$  (Intake) 

1.04 

1.05 

1.02 

1.03 

1.03 

1.04 

4  (exhaust) 

0.95 

1.03 

1.01 

1.00 

1.07 

0.98 

Residual  Fraction^) 

17.3 

15.1 

12.5 

8.5 

11.2 

7.0 

Fuel  Flow  (g/cyc) 

0.0099 

0.033 

0.013 

0.010 

0.034 

0.010 

Fuel  Energy  (J/cyc) 

426 

435 

561 

1428 

444 

1469 

n  (Combustion)  f 

Fuel  Energy  x  n  . 

(J/cyc)canb 

96 

91 

93 

96 

91 

96 

407 

397 

520 

1368 

402 

1406 

Heat  Release 

Q  max  (J/cyc) 

392 

411 

516 

1340 

420 

1440 

Q  max  (f  of  intake) 

Q  max/(FE  x  n _ .  ) 

cano 

91.9 

94.4 

92.0 

94.1 

94.5 

98.0 

96.3 

103.5 

99.2 

98.0 

104.5 

102.4 

e(J/cyc) 

15.3 

-13.9 

4.3 

28.1 

-18.0 

-34.0 

e  norm  {%) 

3.76 

-3.51 

0.83 

2.06 

-4.48 

-2.42 

*  Knock  Limited 


APPENDIX  B 


Charts  and  Graphs 
Ricardo  Engine  Cross-section 
Ricardo  Engine  Cross-section 
Engine  Geometry  and  Friction  Data 
Fuel  Characteristics/Density 
Pressure  Transducer  Calibration 
Fuel  System  Schematic 
Phasing  P-V  Diagram 
PV  Diagram  (scale  and  bias) 

IMEP  vs  Torque  MEP 
IMEP,  Torgue  vs.  m 

A 

(Intake)  versus  t  (exhaust) 
Emissions  versus  $  (exhaust) 


1^ 


HYDRA  MK  IE  ENGINE 
85*7  BORE  X  86  STROKE 

Figure  B-l.  Ricardo  Engine  Cross-seccion 


Engine  Geometry 


Bore 

Stroke 

Compression  ratio 
Correcting  rod  length 


8.57  cm 
8.6  cm 
8.4  :  1 
15.79  cm 


Valve  Opening  and  Closing  Angles 


10 

12° 

BTDC 

1C 

56° 

ABDC 

EXO 

56° 

BBDC 

EXC 

12° 

ATDC 

Friction 


KW 


Rev/sec 

0.9 

1.6 

2.2 

3.2 

27 

w 

57.5 

4.4 

5.8 

6.9 

8.4 

W75 

7T 

85 

95.5 

Figure  B-3.  Engine  Geometric  and  Friction  Data 
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Indolene  Characteristics 


Hydrogen  (J  by  weight) 

12.34 

Carbon  ($  by  weight) 

86.86 

Hydrogen/Carbon  ratio 

based  on  atomic  weight 

1 .69 

Ave.  molecular  weightCC^H^  g^) 

95.83 

Stoichiometric  air/fuel  ratio 

14.37 

Lower  heating  value  BTU/lb. 

18,562.00 

RON 

97.2 

API  gravity 

56.6 

Density  (g/cra^)  of  Indolene  at  Fuel  Temperatures  (°C) 


0.999  x 

CO. 7523  ♦  (15.6 

-  Fuel  Temp(°C)) 

X 

0 

• 

0 

0 

0 

OO 

1—1 

Temp 

Density 

Temp 

Density 

17 

.7504 

24 

.7448 

18 

.7496 

25 

.7439 

19 

.7488 

26 

.7431 

20 

.7480 

27 

.7423 

21 

,7472 

28 

.7415 

22 

.7464 

29 

.7407 

23 

.7456 

30 

.7399 

I 


Figure  B-4.  Indolene  Characteristics 
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Figure  B-6 

Fuel  System  Schematic 


1MEP  gross 


0  20  40  60  80  100  120  140  160 


IMEP(psi)-BMEP (Tor quel +fMEP+PMEP 

Figure  B-9 

Verification  of  engine  1MEF  versus  torque  HEP 
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0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  L.0 

^  A  /  ^  A  WOT 


Comparison  of  IMEP 


gross 


Figure  B-10 

and  BMEP  (from  torque)  versus  air  flow 


» 


Measured  intake  equivalence  ratio 


Mole  fraction  percent  (dry) 
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0  c 
□  c 

A  o 


+  HC  x  10  S 

0  NO  x  lQjr 

x  D 

a 

0  □  \ 

0  \ 

- Tdeal 

CO 

r  * 

□ 

UJ2 

+ 

0 

0  0 

0 

A 

Ideal  02 

+ 

Ideal  CO  / 

exhaust  equivalence  ratio 
0.7  0.8  0.9  1.0  1.1  1.2  1. 

Figure  B-i2 

Exhaust  Emisson  check  versus  ideal  exhaust  emissons 


1500  RPM  0.35  atm 


mm 
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APPENDIX  C 


AFFECTS  OF  VARYING 
MODEL  PARAMETERS 


3.  Crevice  Volume 

4.  Wall  Temperature 

5.  Exponent  (n) 

6.  Intake  Pressure 


Crankang I es 

Figure  C-l 

Affects  of  changing  model  parameter 


1500  RPM  Pj-  0.35  atm 


Crankang les 


Figure  C-2 

Affect  of  changing  model  parameter  C 
1500  RPM  P  »  0.35  atm 


Crankangles 

Figure  C-3 

Affect  of  changing  the  crevice  volume  as  a 
percentage  of  the  clearance  volume 
1500  RPM  P  -  0.35  atm 


Crankangles 

Figure  C-4 

Affect  of  Varying  wall  temperature  (T^) 
1500  RPM  Pj-  0.35  atm 


Loss  [J] 
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Crankang I es 

Figure  C-5 

Affect  of  changing  the  exponent  (n) 
in  the  heat  transfer  equation 
1500  RPM  0.35  atm 


Crankang I es 

Figure  C-6 

Affect  of  varying  intake  pressure  at  1500  RPM 
Data  from  P  *  0.35  atm 
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APPENDIX  D 

•BEST  FIT'  GRAPHS  OF: 

1 )  Combustion  Heat  release 

2)  Cylinder  Energy 

3)  Chemical  Energy  Released 

at  the  six  engine  operating  points  with: 

Ct  -  2.0 

c2  -  0.0 

n  -  0.8 

V  -  1  f  of  V 
or  c 


Notorind  file  J  RICARDO. 24  43  Pal  1  Th]  43 

Firms  file  !  RICARDO.??  43  Cel  1  Id]  43 


ENGINE  RELATED  DATA  J 

FUEL  RELATE!)  DATA  1 

Ce.1  Int.Vlv,Cl.= 

-124,000 

If!  Fuel  T«ee= 

indolenp 

Tdl  Exh.Vlv.Qp.s 

124.000 

Thl  l.HV  CkJ/kd.l3 

0.431E+05 

til  Spark  Loc.= 

-30.000 

Ul  Crev  V.Ccrl= 

0.670E+00 

EXPFRIHFNTAl  CONDITIONS 

• 

• 

INPUT  VALUES  TO  ANAI.YSIS 

4 

» 

Ckl  Inl.P.Catal= 

0.356 

111  Area  at  TDC= 

146.600 

Cal  RPH= 

1000.000 

Cn.l  Exponent3 

0.300 

l.ol  S ui  rl  rate3 

O.OOOE+OO 

I.pI  Itoshni  Cl= 

2.000 

Co.l  Fuel  fd/sl3 

0.824E-01 

Crl  Uoshni  C?= 

0.000 

Csl  Air  Cd/sl= 

0.113E+01 

Ct.l  Resid.frac.= 

17.300 

Cu.7  Ave.  Ft  rind3 

Y 

Cvl  Hall  taeffO 

4 00,000 

C«1  Create  Plot3 

N 

Heat  Release 

A  n  a  1  s  s  i 

s  For  RICARDO  FNGINf 

ANALYSIS  FOR 

AVERAGE 

CYCLE 

Total  Fuel  Fnersu  (J)  J 

476.257 

Total  Charde  at  IVC  (G)  *, 

0.176 

ANALYSIS  OF  IMF  RESUI.TS  FOR  THIS  CYCLE 

Max  in*  pressure!  Paax=10.170  ata  at.  ?3  decrees 

Haxiaua  rate  of  net  heat  release  =  0.123E+02  J/des  at  17  dedrees 

Haxiaua  rate  of  dross  heat  release  =  0.131E+02  J/des  at  17  desrees 

Haxiaua  dross  heat  release  =  0.392E+Q3  Joules 

i.el  91.9  Z  of  total  enerds  input.*  at  102  dedrees 

*  tOZ  of  aax  at  3  dedrees 
t  90Z  of  aax  at  34  dedrees 

thus: 

Initiation  tiae  -  33  dedrees 
and 

Proeadation  tiae  *  31  desrees 


Heat  transfer  at  andle  of  aax  flsross  3  Q.684Et02  Joules 
i.eJ  16.0  Z  of  total  enerdu  input 


Figure  D-lA 
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Motoring  flip  !  RICARDO. 30  43  Cal  t  Col  43 

Finns  file  1  RICARD0.28  43  [r.l  1  CdJ  43 

ENGINE  RELATED  DATA  1  FUEL  RELATED  DATA  : 


Cel  Int.Vlv.Cl.® 

Csl  Exh.Vlv.Qp.® 

Eil  Seark  Loc.® 

CJ.1  Crev  V.Ccc  1= 

-124.000 

124.000 

-45.000 

Q.670EJ0Q 

CM  Fuel  Type® 

Chi  LHV  CkJ/ksl® 

indol ene 
0.43lEt05 

EXPERIMENTAL  CONDITIONS  1 
Ckl  Ini. P. Catal= 

Ca.l  RPM® 

Col  Suirl  rate® 

Cal  Fuel  Cg/sl® 

CsJ  Air  ts/sl= 

Cul  Ave.  Firing® 

Cm3  Create  Plot® 

0.350 

1500.000 

O.OOOE+OO 

0.126E+00 

0.172EL01 

Y 

N 

INPUT  VALUES  TO  ANALYSIS 
Cll  Area  at  TDC® 

Cnl  Exponent® 

I>1  Uoshni  Cl® 

Crl  Uoshni  C2= 

Ctl  Resid.frac.® 

CvJ  Uall  tapCK.l* 

« 

a 

146.600 

0.800 

2.000 

0.000 

15.100 

400.000 

Heat  Release 

A  n  a  1  s  s 

is  For  RICARDO  ENGINE 

ANALYSIS  FOR 

AVERAGE  CYCLE 

Total  Fuel  Energy  (J)  1 

435.224 

Total  Charge  at  I  VC  (G>  1 

0.174 

ANALYSIS  OF  THE  RESULTS  FOR  THIS  CYCLE 

Maxima  pressure!  Paax=12»847  ate  at  16  degrees 

Maxima  rate  of  net  heat  release  *  0.119EE02  J/deg  at  9  degrees 

Maxima  rate  of  dross  heat  release  =  0.130E+02  J/deS  at  9  degrees 

Maxima  gross  heat  release  *  0.411E+03  Joules 

i.el  94.4  Z  of  total  energy  ineuti  at  114  degrees 

t  10Z  of  aax  at  -8  degrees 
t  90Z  of  aax  at  25  degrees 

THUS! 

Initiation  tiae  =  37  degrees 
and 

Propagation  tiae  =  33  degrees 


Heat  transfer  at  angle  of  aax  Osross  ®  0.790E+02  Joules 
i .el  18.1  Z  of  total  energy  input 

Figure  D-2A 
1500  RPM  P  -  0.35  atm 


s  [J] 
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Hntorind  file  :  RICAR00.45  43  Tel  t  Chi  43 

Firind  file  1  RICARD0.43  43  Cel  1  Cdl  43 


ENGINE  RELATED  DArA  l 

FUEL  RELATED  DATA  : 

Cel  Int.Vlv.Cl.® 

-174.000 

CM  Fuel  Twee® 

indolene 

Cdl  Exh.Vlv.OP.® 

124.000 

Chi  LHV  CkJ/kdl® 

0.43!Et03 

Ci 1  Seark  Lor.= 

-40.000 

f  J1  Crev  V.Cccl® 

0.470E+00 

EXPERIMENTAL  CONDITIONS 

« 

• 

INPUT  VALUES  TO  ANAl.YSTS 

• 

♦ 

Ckl  Inl.P.Cat»T= 

0.434 

Cl!  Area  at  TDC® 

144.400 

Cal  RPN® 

1500.000 

Cnl  Exponent- 

0.800 

Col  Suirl  rate® 

0.000E+00 

Cel  Hoshni  C.t= 

2.000 

Cal  Fuel  Cd/sl® 

0.143E+00 

Crl  Hoshni  C7= 

0.000 

Csl  Air  Cd/sl® 

0.230E+01 

Ctl  Resid.fr»c.= 

12.500 

Cul  Ave.  Firind® 

Y 

Cv!  Hall  tapfKl® 

400.000 

Cm!  Create  Plot® 

N 

Heat  Release 

Anal  h  s  i 

s  For  RICARDO  ENGINE 

ANALYSIS  FOR 

AVERAGE 

CYCLE 

Total  Fuel  Enerdu  ( J>  : 

541.101 

Total  Charse  at  IVC  (G)  l 

0.775 

ANALYSIS  (IF  THE  RESULTS  FOR  THIS  CYCLE 

Maxima  eressurel  p»ax=18.825  ate  at  13  decrees 

Maxima  pate  of  net  heat  release  =  0.190E402  J/ded  at  5  decrees 

Maxima  rate  of  dross  heat  release  =  0.208ET02  J/deS  at  5  decrees 

Maxima  dross  heat  release  =  0.514E+03  Joules 

i.el  92.0  X  of  total  enerds  inputi  at  113  dedrees 

t  10Z  of  tax  at  -10  desrees 
t  90Z  of  aax  at  14  dedrees 

thus: 

Initiation  tiee  =  30  desrees 
and 

Propagation  tiae  =  24  decrees 


Heat  transfer  at  andle  of  aax  Qdross  =  0.999E+02  Joules 
i.el  17.8  Z  of  total  enerdx  input 

Figure  D-3A 
1500  RPM  P_j.-0.44  atm 


CJ] 


Crankangles 

Figure  D-3B 
1500  RPM  Px-  0.44  atm 
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Hotorins  file  i  RICARBU.8 

43  Cal  1 

Ch;  43 

Firing  file  2  RICARDO.l 

43  Ed  1 

Cril  43 

ENGINE  RELATED  DATA  i 

FUEL  RELATED  DArA  ! 

Cel  Int.Vlv.Cl.® 

-124.000 

Cf.T  Fuel  Twee* 

indoiene 

CS.T  Exh.Vlv.0p.* 

124.000 

Ch.T  I.HV  CkJ/ksl* 

0.431E405 

Ci.T  Seark  Loc.= 

-22.000 

Evil  Crev  V.Cccl* 

0.470E+Q0 

EXPERIMENTAL  CONDITIONS  J 

INPUT  VALUES  TO  ANALYSIS 

• 

• 

CkJ  Ini .P. fatal* 

1.008 

Ill  Area  at  TBC* 

146.6CG 

Cal  RPM* 

1500.000 

Cn]  Exponent* 

0.800 

Col  Suirl  rate= 

0.000F+00 

Cel  Uoshni  Cl* 

2.000 

Cal  Fuel  Cs/sl* 

0.414E+00 

Crl  Uoshni  C2* 

0.000 

Csl  Air  Cs/sl* 

0.579E+01 

Ct.l  Resid.frar.* 

8.500 

Cul  Ave.  Firins* 

Y 

Cvl  Wall  tapCKl* 

400.300 

Cul  Create  Plot* 

Heat  Release 

A  n  a  1  n  s 

is  For  RICARDO  FNGINE 

ANALYSIS  FOR 

AVERAGE  CYCLE 

Total  Fuel  EnerSs  <J)  2 

1428.100 

Total  Charse  at  JVC  (G)  ! 

0.542 

ANALYSIS  OF  THE  RESULTS  FOR  THIS  CYCLE 

Maxiaua  pressure!  Paax=42.547  ata  at  18  desrees 

Maxiaua  rata  of  not  heat  release  *  0.551E+02  J/des  at  10  desrees 

Maxiaua  rate  of  dross  heat  release  =  0.590E+02  J/des  at  10  desrees 

Haxiaua  Sross  heat  release  =  0.134E+04  Joules 

i.el  94.1  Z  of  total  enerss  ineut»  at  100  desrees 

*  102  of  oax  at  -3  desrees 
<  902  of  aax  at  20  desrees 

thus: 

Initiation  tiae  *  19  desrees 
and 

ProeaSation  tiae  *  23  desrees 


Heat  transfer  at  ansle  of  aax  Osross  *  0.187E+03  Joules 
i.et  13.1  2  of  t-ot-sl  enerss  ineut 

Figure  D-4A 
1500  RPM  WOT 


A 
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Crankang les 

Figure  D-4B 
1500  RPM  WOT 
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Motoriins  file  *  RICARDO. 40 

43  Cal  ? 

Chi  43 

Firms  file  :  RICARDO, 38 

43  CM  i 

Cdl  43 

ENGINE  RELATED  DATA  : 

FULL  RELAYED  DAfA  '. 

Cel  Int.Vlv.Cl .® 

-124.000 

Cf3  Fuel  Tspe® 

indolene 

Cd.T  Exh.Vl v.Oe.® 

124.000 

Chi  LHV  CkJ/kdl® 

0.4312+05 

CM  Seark  Lnc,® 

-45.000 

CJ1  Crev  V.CccI® 

0.470E+00 

EXPERIHENTAL  CONDITIONS  i 

INPUT  VALUES  TO  ANALYSIS 

« 

Ckl  Inl.P.Catdl® 

0.350 

Cl.l  Area  at  TDC® 

144.400 

Cal  RPN® 

2500.000 

CnI  Exponent® 

0.800 

Col  Swirl  rate® 

0.000E+00 

Cp.I  Uoshni  Cl® 

2.000 

Col  Fuel  td/sl* 

0.215E+00 

Cr.l  Wnshni  C2= 

0.000 

Cs.l  Air  Cd/sl® 

0.297E+Q1 

Ctl  Resin, frac.® 

11.2C0 

Cul  Ave.  Firind® 

Y 

Cvl  Wall  teeCIQ® 

400.000 

Cwl  Create  Plot= 

N 

Heat  Releas 

e  Analwsi 

s  For  RICARDO  ENGINE 

ANALYSIS  F  0 

R  AVERAGE 

CYCLE 

Total  Fuel  Enerdw  (J) 

« 

:  444.052 

Tota2  Charse  at  JVC  (fi)  J 

0.172 

ANALYSIS  OF  THE  RESULTS  FOR  THIS  CYCLE 

Maxieue  pressure!  Piax=J3.231  ate  at  IS  decrees 

Haxj»»  rate  of  net.  heat  release  =  0. 1242+02  i/ded  at  10  decrees 

Haxiaua  rate  of  dross  heat  release  =  0.134E+02  J/ded  at  10  dedrees 

Haxiaue  dross  heat  release  *  0.420E+03  Joules 

i.e.  94.5  Z  of  total  enerds  ineuti  at  115  dedrees 

t  10Z  of  eax  at  -tl  decrees 
*  90Z  of  >ax  at  24  decrees 


thus: 


Initiation  tide  ®  54  dedrees 
and 

Propagation  tite  -  35  desrees 


Heat  transfer  at  ansle  of  tax  Odross  -  0.752E+02  Joules 
i.e:  14. 9  Z  of  total  enerds  input 

Figure  D-5A 


2500  RFM  P^O.35  atm 


ss  [J] 
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-150  -100  -50  0  50  100  150 

Crankangles 


Figure  D-5B 
2500  RPM  Pr»  0.35  atm 
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Motorfod  fiJ#  :  ricardo.18  43  tal  )  tbl  43 

Firing  file  >  ricardo.15  43  tel  t  Edl  43 


ENGINE  RELATED  DATA  5 

FUEL  RELATED  DATA  1 

tel  Int.Vlv.CU> 

-124.000 

tf 1  Fuel  Tvre> 

indotene 

CsJ  Exh.Vlv.0e.> 

124.000 

thl  LHV  CkJ/kd> 

0.431E405 

til  Searit  loo 

•30.000 

fJl  Crev  V.Cccl* 

0.470E+ 00 

EXPERIMENTAL  CONDI  HONS 

1 

INPUT  VALUES  TO  ANALYSIS  1 

Tkl  Inl.P.tata> 

1.003 

til  Area  at  TDC* 

140.400 

tel  RPH> 

2500.000 

tnl  Exponent- 

0.800 

to]  Swirl  rate* 

O.OOOEtOO 

EpI  Noshni  Cl> 

2.000 

Eel  Fuel  ts/sl» 

0.710E+00 

trl  Noshni  C2> 

0.000 

tsl  Air  Cd/sl* 

0.993E+01 

Ctl  Resid.ffac.a 

7.000 

tul  Ave.  Firins* 

V 

tvl  Hall  tertJCl- 

400.000 

twJ  Create  Plot* 

n 

V;  . 

Heat  Release 

A  n  a  1  v  s  i 

s  For  RICARDO  EN8INE 

V 

ANALYSIS  FOR 

AVER  A8E 

CYCLE 

Total  Fuel  Enordu  <J)  1 

1409.13? 

Total  Charde  at  IVC  (6)  J 

0.549 

ANALYSTS  OK  THE  RESULTS  FOR  THIS  CYCLE 

tori pressure!  P»ax>4?.45A  at*  at  18  dedrees 

Maximo  rata  of  nat  heat  rtliiM  ■  0.574E+02  J/ded  at  tl  dedrees 

tori  mm  rata  of  dross  heat  release  *  0.411E4O2  J/ded  at  11  degrees 

Maxieuo  dross  heat  release  *  0.144E+04  Joules 

i 98.0  Z  of  total  enerdv  input.  at  97  desrees 

t  10Z  of  eax  at  -2  desrees 
t  90Z  of  tax  at  23  dedrees 

must 

Initiation  tine  *  28  dedrees 
and 

Propadation  tiae  »  25  dedrees 


Meet  transfer  at  andle  of  aax  Qsross  •  0.173E+03  Joules 
i.el  11.8  Z  of  total  enerds  ineut 


Figure  D-6A 
2500  RPM  WOT 


Least  Squares  Error 
Calculations  to  Determine 
the  Model  Parameters 
Yeilding  "Best  Fit" 
Heat  Release 
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APPENDIX  F 


GRAPHS  FOR  COMPARISONS 
OF  LOADS  AND  SPEEDS 

1 )  Cylinder  Energy 

2)  Chemical  Energy  Released 

3)  Combustion  Heat  Release  Profile 

with  Model  Parameters: 

C1  -  2.0 

Cg  «  0.0 

n  -  0.8 

V  •  1  H  of  Vn 
cr  o 


s  CJ] 


Crankangles 

Figure  F-3 

of  load  at  1500  RPM  with  Pj-  0.35, 


-150  -100  -50  0  50  100  150 


Crankangles 

Figure  F-4 

Comparison  of  load  at  2500  RPM  with  P^*.35  and  WOT 
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Appendix  G 

Fuel  -  Air  Flow  Calculations 

1 .  Calculations 

2.  Air  Flow  Temperature  Corrections 

3.  Fuel  Flow  Chart 


APPENDIX  G 


Calculations  of  Air  and  Fuel  Flow 


Air  Flow  Rate 

The  scale  depends  on  the  maximum  position, 

-c 

(1)  For  top  inclined  position  -  7.50  x  10 

-5 

(2)  For  second  position  -  3.00  x  10 

o 

Air  flow  (m  /s)  -  maximum  reading  x  scale  factor  x  temperature 
correct 1  cxi  (figure  G-2) 

*lr  flan  (f)-  air  n»  (f)  x  (  —^1—^  )  x  <^>1.293  x  103 

where  P(LFE)  •  Barameter  reading  -  pressure  drop  across  the  viscous  flow 
meter  laminar  element 


Fuel  Flow  Rates 

Fuel  Flow  (g/s)  -  measured  volune  flow(cm3/mln)  x  density  of  Indolene(flg.B-4) 

60  x  103 


Equivalence  Ratio  Intake 


Fuel  Flow(g/s)  1  /F. 

Air  Flow  g/sec  x  VTff  V 


Equivalence  Ratio 


Fuel  Setting  80  /  100 


Figure  G-3A 
Fuel  flow  rate 
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